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ABSTRACT 


-  A  statistical  measure  of  the  performance  of  an  adaptive  beam-forming  technique  for  a 
multisensor  linear  array  has  been  developed.  The  sample  complex  vector  signal  process,  Z(/*) 
as  observed  at  sample  times  (/*}  €  {0,7}  is  used  to  calculate  a  positive,  definite,  maximum 
likelihood  estimate  (M.L.E.)  /{'of  the  signal  process  covariance  matrix  R.  Then,  an  estimated 
measure  of  the  performance  of  the  array  in  terms  of  its  ability  to  suppress  an  interfering  signal 
is  developed  from  an  a  priori  specified  process  matrix  E(K),  its  sample  M.L.E.  estimate  /{'and 

a  form  of  the  Wishart  distribution.  In  particular,  such  an  estimated  measure  of  performance  in 
terms  of  suppressing  an  interfering  signal  with  an  angle  of  arrival  ^OA)  of  0j,relative  to  a  beam 
AOA  of  0,  is  the  statistic  represented  by  the  absolute  difference  of  the  Lj'norms*  { I  H/r’Uj* 

-  Iln'lljl }  •  The  estimate  /r1  =  V*  (0,)  R'  K(0,)  ,  it's  exact  value  is  /r‘  =  F*(0,)  R  '  K(0,)  and 
R)R  arc  functions  of  K(03)  and  V*  (03). 
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I.  BACKGROUND. 

The  research  undertaken  here  is  a  statistical  study  of  the  performance  properties  of  an 
adaptive  beam-forming  algorithm  when  utilized  to  process  stochastic  signals  as  received  by  a 
multi-element  linear  array. 

A  performance  measure  is  defined  which  represents  the  degree  to  which  the  adaption  process 
minimizes  off-boresight  signal  interference  from  a  single  source. 

It  is  well  known,  particularly  in  acoustics,  that  multipath  signal  interference  in  a  dispersive 
medium  such  as  salt  water  presents  many  difficulties  to  the  detection  and  acquisition  of  wanted 
signals  when  in  the  presence  of  unwanted  or  interfering  signals.  In  fact,  even  in  the  case  of  a 
single  interferer,  it  is  often  very  difficult  to  supress  such  signals  at  low  signal-to-noise  (SNR) 
ratios  without  the  processing  enhancement  as  offered  by  adaptive  beam-forming  techniques. 

Most  of  the  existing  adaptive  signal  processing  techniques  deal  with  the  deterministic, 
expected  or  steady  state  properties  of  stochastic  signal  processes  which  enable  one  to  defer  or 
even  ignore  the  distributional  aspects  of  the  signal  processes  being  considered.  Therefore,  since 
most  interest  and  concern  has  been  given  heretofore  to  beam-forming  techniques  using  assumed 
expected  values  of  signal  processes  parameter  estimates  rather  than  to  the  distributional  aspects 
thereof,  it  was  decided  that  improvements  in  the  analysis  of  beam-forming  techniques  could  best 
be  achieved  by  a  rigorous  statistical  analysis  of  the  problem.  Thus,  major  emphasis  was  placed 
on  development  of  the  distribution  of  the  estimated  measure  of  supression  of  off  boresight 
signals. 

It  has  long  been  recognized  that  the  distributional  aspects  of  statistical  sampling  theory' 
could  be  adapted  to  provide  quantitative  measures  of  performance  for  physical  processes  such 
as  the  beam  forming  of  signals  from  an  array  of  sensors.  In  some  cases,  such  as  for  multi-var¬ 
iate,  (vector)  Gaussian  stochastic  processes,  many  of  the  essential  sampling  distributional 
properties  and  characteristics  were  developed  and  published  by  John  Wishart3  of  Rothamsted 
Experimental  station  in  England,  circa  1928.  Then  in  those  lean,  worldwide  depression  years 
of  1932  and  1933,  Wishart  and  Bartlett3-4 ,  while  at  Cambridge  University,  utilized  a  character¬ 
istic  function  approach,  made  possible  by  the  work  of  Ingram5 ,  to  rigorously  established  the 
existence  of  a  class  of  distributions  which  became  known  later  as  Wishart  distributions.  In 
1938,  the  noted  Chinese  mathematician,  P.  L.  Hsu* ,  generated  an  elegant  proof  for  the  exist¬ 
ence  of  a  general  form  of  the  Wishart  distribution  by  mathematical  induction. 

In  1956  Cramer1  and  in  1963,  Fisz*  exemplified  this  form  of  the  Wishart  distribution,  citing 
the  earlier  works  of  Wishart  and  Bartlett  and  the  existence  proof  of  Hsu,  respectively.  With  the 
advent  oflater  works  in  the  vector  calculus,  linear  algebra,  matrix  theory  and  modern  abstract 
analysis,  multivariate  statistical  analysis  became  a  very  powerful  tool  in  the  mathematical 
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modeling  of  multidimensional  (vector)  stochastic  physical  processes  such  as  encountered  in  the 
adaptive  beam  forming  of  passive  signals,  received  from  an  array  of  sensors  consisting  of  a  large 
number  of  elements. 

In  1984,  Anderson9  of  Stanford,  published  the  second  edition  of  his  "Encyclopedic”  text 
on  multi-variate  analysis  in  which  he  rigorously  developed  all  of  the  properties  of  the  several 
classes  of  Wishart  distributions  both  in  the  real  and  complex  domains.  Thus,  it  is  the  Wishart 
distribution  and  the  various  transformed  versions  thereof  which  provide  the  statistical  founda¬ 
tions  for  the  work  to  be  presented  in  this  report. 
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2.  INTRODUCTION. 


In  the  research  results  presented  here,  a  statistical  measure  of  performance  for  an  adaptive 
beam-forming  technique  and  its  associated  distribution  function  has  been  developed  for  the  case 
of  a  complex  Gaussian  signal  process.  The  statistical  measure  is  a  sample  estimate  of  the  degree 
to  which  and  adaptive  beam-forming  process  fot  *  linear  array  can  minimize  the  ofT-axis  inter¬ 
ference  of  a  plane  wave  propagated  from  a  single  interferer.  Most  of  the  existing  signal 
processing  references  are  concerned  only  with  the  expected  value  of  requisite  statistical  estima¬ 
tors.  That  is,  the  existing  work  has  assumed  only  that  the  sample  estimators  of  interest  exist 
and  are  unbiased.  What  is  not  well  understood  by  many  signal  processing  practitioners, 
however,  is  how  system  parameters  and  the  number  of  sensors  effect  the  distribution  of  the 
estimators.  When  complete,  this  research  can  be  utilized  to  provide  answers  to  the  following 
types  of  questions: 

•  Given  a  noise  model,  interference  mode!  and  sensor  array  size,  what  is  the  optimal  sampling 
method  in  order  to  obtain  a  specified  probability  of  detection  while  at  the  same  time 
utilizing  a  minimum  amount  of  computing  resources? 

•  What  arc  the  tradeoffs  being  made  when  a  system  design  selects  block  averaging  over 
exponential  averaging ?  Will  system  performance  improve  and  if  so  by  how  much? 

•  It  is  well  known  that  the  noise  sampled  at  adjacent  sensors  is  not  statistically  independent. 
What  are  the  effects  of  this  anisotropy  on  the  distribution  of  the  beam-former  output  as  a 
function  of  sampling  technique? 

•  For  arrays  with  a  very  large  number  of  sensors  (  >  10’)  a  long  sampling  time  would  be 
required  to  form  a  good  estimate  of  the  covariance  matrix.  In  this  case,  how  docs  one 
develop  a  strategy  to  combine  sensors  to  produce  the  best  improvement  in  the  estimation 
of  the  covariance  matrix  while  keeping  within  the  desired  integration  times? 


The  distribution  of  the  performance  estimate  (in  terms  of  the  inverse  covariance  matrix 
estimate)  enables  the  specification  of  confidence  bounds  for  the  process  performance  measure 
as  a  function  of  the  number  of  samples  and  an  a  priori  selected  level  of  significance.  These 
results  arc  valid  for  any  signal  processing  system  operating  in  a  low  SNR  environment  and 
where  statements  of  precision  for  requisite  estimates  arc  most  desirable. 
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3.  THE  SIGNAL  SAMPLE  SPACE. 

Consider  an  n  element  ( n  ^  2),  multi-sensor  linear  array  with  a  spacing  between  sensors  of 
one  half  the  wave  length  (  ^/2  ) .  The  sample  signal  process  Z(/*)  e  O',  k  =  1,2, ... ,  M  is 
considered  to  be  a  complex,  non-zero  mean,  wide-sense  stationary,  Gaussian  vector  process 
with  sample  components,  /,(/*),  ... ,  z„(/*),  for  each  {/*}  e  [0,7],  T<  oo,  k  =  1,2, ... ,  M. 

Thus,  the  sample  vector 


and  Z*  (/*)  =  [zj(/*)  22 (/*) ...  zn(4)] 


Then,  for  each  /,  e  [0,7];  k  =  1,2, ... ,  M  ,  the  complex  sample  vectors  Z(/,),  Z(/3), ,  Z(iM)  arc 
assumed  to  be  independent  identically  distributed  (I.I.D.)  over  a  temporal  sample  interval 
[0,7]  for  7'<  oo.  Perhaps  a  less  restrictive  condition  than  general  independence  is  to  consider 
the  sample  vectors  as  random  variables  from  a  quasi-stationary,  Gaussian  sample  process  which 
is  (strict  sense)  stationary  over  the  finite  sample  interval  [0,7],  but  temporally  uncorrelated 
through  at  lease  the  second  degree.  In  general,  a  second  degree,  uncorrelated  sample  process 
is  a  process  in  which  the  sample  values  about  the  mean  are  uncorrclated  such  that 

£[(4^)%,  =  0for(M=  1.2 . M. 

A  principle  reason  for  meeting  such  a  sampling  condition  is  that  the  sample  estimates  of  the 
process  parameters  which  define  the  covariance  matrix  are  statistically  consistent  and  hence  so 

A 

is  the  covariance  matrix  R. 


•  •  *  •  *  V  V  ».  \  o 
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4.  COMPLEX  MULTIVARIATE  GAUSSIAN  CASE. 


Consider  a  sample  set  of  size  M  of  second  degree,  uncorrelated,  Gaussian  random  vectors 
{Z,  ,  Zj ... ,  ZMj  e  O  where  Z ,  =  Z(t ,),  i  e  [1,A /],  as  previously  defined.  Then  M  (n  element) 
I.I.D.  complex  vectors  have  the  joint  probability  density  function, 


/(Z„  Z, . Zw)  =  f|(,TL)|/!|-''i 


*■1 

7rr-\R\^f  st-m-*?*-' 


where  the  complex  process  means  0,  =  E(Z)  and  the  covariance  matrix  R  e  M„(n  x  n  complex 
matrix)  is  a  nonsingular  llermitian  matrix  such  that  {Z  —  ©,)*  R(Z  -  0,)  >  0  for  all  non-zero 
(Z  -  0.)  e  O,  and  hence  R  is  positive  definite.  Therefore,  the  matrix  R  1  e  Mn  exists  and  is  also 
a  positive  llermitian  matrix. 
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5.  MAXIMUM  LIKELIHOOD  ESTIMATES  OF  0  R. 

It  can  be  demonstrated  that  if  2,2,, ... ,  2*  constitutes  a  sample  of  size  M  of  an  n  element 
1.1. D.  complex  Gaussian  vector  process  where  each  complex  vector  is  N(Q,  R )  with  M>  n  , 
then  the  maximum  likelihood  estimators  of  ©  and  R  are: 


m  m 

A  1  V1  A  I  A  A 

©,=-77  \zkR  ~  ~  }jzk  -  e,)(zk  -  o2y 


For  computational  purposes  we  write 


Tzik 

km  |  km  1 

(5.2) 

M 

\  \'  •  A  A  * 

-jjlZkZk-M®2®2 

Ar«! 

(5.3) 

Z{ik)  =  Zk  =  and  7*{ik)  =  7k  =  [?,*  z2k  2nk] 


M  M 

(5-5) 

km  1  *-l 

M  M 

*-l  *-1 

(5.6) 

and  hence,  the  sample  covariance  matrix  is 


M  M 


R,=  [ 7,1  ;  4 -  -Jr  Iz*z;* -  Zz'*Iz/- 


*-i  *-i 
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6.  THE  WISHART  DISTRIBUTION. 


In  1928  John  Wishart  prublishcd  the  first  of  several  papers 1X*  in  multivariate^statistical 
analysis.  Anderson*  shows  that  if  the  sample  covariance  matrix,  R  —  (M  -1)5  =  £  (Xt  —  A')', 
where  X„ .... ,  X^M>n)  are  n  component,  real  valued,  statistically  independent  vectors,  each 

A 

with  distribution  N(n,  R)  then  R  has  the  density 


g(R\R,M,n) 


|  ^  |  )l2(M-n-2)g-\l2tr(R'lR) 

I'M"-" -5^- 1 R |  '“"-"fjrC  1  /2(M  -0] 


(6.1) 


for  R  positive  definite  and  0  otherwise.  The  expression  (5.1)  is  the  original  Wishart  result  in  real 
vector  form.  For  Z,,  Z}, .... ,  ZN  normally  distributed  complex  vectors,  each  with  mean  O  and 

A 

Hermitian  covariance  matrix,  R,  the  sample  estimate  R  of  the  covariance  matrix  R  has  the 
complex  Wishart  density 

g(R  e  Mn)  |  RM,n)  = - LJ - T -  (6.2) 

lm  1 


where  R  =  £(Z»  -  0,)(Z1  -  0Z)*,  0,  =  ~  J]Z*  ,  M=  N  -1  ,  and  (R  and  R)  e  M„  are  Hermi- 

1-1  _  M.  *>i 

tian  positive  definite  matrices,  since  Z*  RZ  >  0  and  Z*  RZ  >  0  for  all  non-zero  ZeO.  The 
matrices  R  and  R  have  the  additional  property,  that  if  the  number  of  sample  vectors  M  arc 
greater  than  n  the  number  of  components  of  each  vector,  then  the  probability  is  unity  of 

A  A 

drawing  a  sample  so  that  R  is  positive  definite  (Z*  /?Z)  >  0  and  on  the  other  hand,  if  M  <  rt 
(i.e.,  the  number  of  sample  vectors  is  less  than  the  number  n  of  the  components  of  which  each 

A 

vector  is  composed)  then  it  can  be  demonstrated  that  R  does  not  have  a  density,  but  is  does 
have  a  well  defined  cumulative  distribution  function  (c.d.f.)  (monotone  non-decreasing). 


G(R  |  R,  M,  n)  = 


0  \R\<,0 

0Za<,  1  |£|  >0 


(6.3) 


A  A 

For  R  and  R  Hermitian  positive  definite  it  follows  that  R~'  and  R-' 
positive  definite.  In  addition  by16 


are  also  Hermitian 


Theorem  1 

If  //  is  a  Hermitian  matrix,  then  x*11x  is  real  for  all  xe  O. 
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Proof: 


Write  (x*Ilx)  =  (x *Hx)*  =  (*)*(//)*( x*)*  =  x'll+x)  =  (x'llx),  i.e., 
conjugate  and  is  therefore  real  by  definition. 


(x*Hx)  equals  its  complex 


6-2 
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7.  BEAM-FORMING  MEASURE  OF  PERFORMANCE. 

In  adaptive  beam-forming  problems,  an  appropriate  measure  of  performance  of  the  beam- 
froming  technique  is  the  degree  to  which  the  beam-former  output  approaches  an  ideal  supres- 
sion  response  to  unwanted  signals  in  an  otherwise  “noise  free”  environment.  By  noise  free  it  is 
meant  in  the  absence  of  external  and/or  internal  system  noise,  deterministic  or  random. 

Thus,  an  ideal  measure  of  performance  in  the  sense  of  suppressing  the  unwanted  signals  of 
a  single  interferer  at  an  AOA  of  0,  by  the  beam-forming  processing  at  an  AOA  of  0,  ,  is 

p~'  «  V*(0\)R~\0i)V(0\)  (7-05) 

Then  the  sample  estimator  p of  p *•  and  its  destributional  properties  with  resepct  to  p  ' 
are  defined  as  follows. 

Consider  an  interfering  signal  with  an  AOA  of  03  relative  to  a  linear  array  with  a  beam 
formed  at  an  AOA  of  0,.  For  /?*'(0 ,)—  [r,/0j)3  Hermitran  and  positive  definite,  it  follows 
Pi}(0„  9j)  =  y*(9,)R  '(9,)V(9I)  >  0  (postivc  definite)  and  Teal  by  Theorem  1.  F*(0),  V(9)  are  the 
usual  phase  vectors  used  in  conventional  beam  forming.  Because  pi,'(0„  0})  is  a  function  of  the 
number  of  time  samples  M  for  any  AOA's  0„03  it  can  be  demonstrated  that 

(7» 

and 

M-oo  =  0  (7-2) 

Hence,  p*1  converges  in  expectation  to  p  '  and  is  also  a  consistent  estimate  by  (7.2). 

For  notational  sake,  we  will  use  (%„)  —  R-]  and  (if,,)  =  R'K  The  sample  statistic  p-M'  can  be 
written  as  a  matrix 


£,'-£[  f7-3' 

U 

which  is  positive  definite  and  thus  real  by  Theorem  I. 

A 

Deriving  the  distribution  of  R  l  from  the  complex  Wishart  density  given  by  (5.2)  and 

A  A 

computing  the  Jacobian  of  the  transformation  (/?  =  R  '),  which  is  the  determinant  of  the 
rt3  x  n3  positive  definite  matrix 
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m-= 


Sr!(j 


bn kj 


=  \R  I 


—In 


(7.4) 


Consequently, 

g(B)  =  |  (R,M,n)  I  *  I"2"  =  g(R']  |  R~] ,M,n ) 

A 

if  R-'  is  positive  definite  and  0  otherwise. 


(7.5) 


The  density  distribution  for  the  sample  statistic,  p~'  can  then  be  obtained  from  (7.5)  for 
B  =  p  x  =  V*  (0,)K-'K(0,)  and  calculating  J{B )  =  ./[K^fl-'F^,)]  =  J[bJ  where  b„ 

-  [»?,,v,(0|)v;(0,)]  .  These  calculations  yield 


A  *  A 


g(b  |  r,M,")J(B)  =  g(V  (0,)K  V{0\)  I  PM,n) 
for  p  1  ;>  0  and  0  otherwise. 

Using  the  Wishart  density  function  (7.6)  it  is  found  that  the  sample  statistic, 


(7.6) 


n 

£ltyl2vMW) 


(7.7) 


with 


n 


(7.8) 


has  density 

g(b  |  P,M,n)J(B) 
then  it  follows  quickly  that, 


H 


sir'll: I \\p-'h,M,n)  for  llp-'ll^O] 


0  otherwise 


(7.9) 


P[\  Rp-Ilh  -  !IP-Ilh!  ^c(A/,  a)]  =  fW(||p-'ll2  \  fip''ll2.w.J=  l  -o 

Jo 


(7.10) 


or 


ftjvW.  «) 

r{[fp-'h  -  llP_’lb]  ^  « \M,  a))  =  dG(\\p-%\  lip"'1 

Jo 


112.  M.n 


)  =  1  -a 


(7.11) 


7-2 


*>-*-•  mJr  V. 


.*>_  *.->  V  *»i  >*-»  *J*  ' 


3 


£ 


X 

«/\ 


y. 


*  * 

>1 

•  A 

V' 


I 

is 

•S 

<5 


V, 

V 

*  *V>  “_A  •_*' 


Where  M  is  the  number  of  time  samples  and  0  <  a  <  0.10  is  an  a  priori  selected  level  of  signif¬ 
icance.  The  cumulative  distribution  function  represented  by  (7.1 1)  yields  the  probability 
1  -  a,  that  H/r'Ili  will  fall  into  ar  interval  (confidence  interval)  of  width  2  c(a,  M)  with  midpoint 
I  P-'II.- 
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8.  A  QUICK  LOWER  BOUND  IN  PROBABILITY. 

Many  times  a  useful  lower  bound  in  probability  can  be  obtained  by  the  use  of  Chebyshev's 
Inequality  Consider  the  sample  statistic  jjp-'ljj  with  <  oo  and  <Ta(IIP_1lla)  <  oo  as  before. 
Then  by  Chebyshev's  inequality  and  the  sample  estimate  d3M, 

P([!!p_,il2  -  £1Ip_1II2]2  £  cJ)  -  1  -  (8.1) 

This  bound  can  be  quite  sharp  tn  this  case,  since  ®UIIp_,Hi)  ”  Another  very  useful 

inequality  for  a  quick  estimate  of  the  statistical  magnitude  of  the  random  variable  ||p-'|i,  is  given 
in  terms  of  the  Maikov  inequality 

/>[llp”,|li  ^  a]  =  (8-2) 

where  a  >  0  and  ||p-'||,  £  0,  regardless  of  the  actual  distribution  of||p  '!!„  which  of  course  has 
been  demonstrated  to  be  Wishart  by  expressions  (7.9)  and  (7.10),  respectively.  For  large 
numbers  of  temporal  samples  M  <  oo,  both  (8.1)  and  (8.2)  can  often  give  very  useful  results 
quite  rapidly.  The  Markov  inequality  has  the  additional  advantage  that  if  the  sample  vector 

A 

components  used  to  calculate  R  1  are  temporally  correlated  as  well  as  spatially  correlated  across 
an  array,  Fllp  'llj  <  oo  can  be  demonstrated  to  be  0(1/Af)  and  is  thus  asymptotically  unbiassed 
in  expectation. 

If  it  can  be  demonstrated  IJ  that  the  estimated  variances  and  corvariances  have  temporal 
component  correlations  of  less  than  for  the  number  of  temporal  samples  M ,  where 

l<Af<oo  for  each  of  the  sensors,  2  <,  n<  oo,  in  the  array,  then  the  Lj  sample  estimates 
such  as  ||pM'||j  are  \\pjj  -  p  %  arc  stoclastically  convergent,  and  asymptotically  unbiassed. 

0(1 /A  f). 
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9.  CONCLUSIONS. 


•  Optima!  sampling  occurs  in  the  adaptive  beam-forming  estimation  process  when  the 
number  of  temporal  samples  M  ^  n,  the  number  sensors,  2.  If  M  <  n,  see  e  g.  (6.3),  the 
probability  density  in  the  form  of  (7.9)  does  not  exist. 

A  \ 

•  Block  averaging  as  required  in  the  sample  estimation  of  R-'  and  hence  pjJ  provides  statis¬ 
tically  consistent  estimates  of  the  process  parameters  necessary  to  the  efficient  interference 
supression  of  unwanted  signals  by  the  beamformer  process. 

•  Exponential  averaging  docs  not  provide  numerically  or  statistically  consistent  esti¬ 
mates  of  the  requisite  process  parameters  either  for  deterministic  or  stochastic  sample 
process. 

•  Spatially  correlated  signal  process  samples  between  adjacent  sensor  elements  of  an  array 

*  A 

arc  accounted  for  in  the  estimation  of  and  hence  the  subsequent  density  distrib¬ 

ution  given  by  7.9. 

•  '  In  the  case  of  a  large  number  of  sensors  n,  (n  >  10')  and  the  requirement  that  the  number 

of  temporal  samples  M^n,  for  the  obtainment  of  statistically  consistent  supression 
performance  estimates,  integration  times  consistent  with  specific  improvements  in  SNR 

ratios  should  pose  no  unusual  requirements  for  computational  agility.  V-c>, ,,,  (s  • 

'  •'  -  r  •  •  <■.  i  .  «r  ■ 

•  A  derived  form  of  the  complex  Wishart  distribution  for  the  sample  estimator  of  an  array's 
performance,  G{||/>  '(!,} ,  provides  a  probability  measure  for  the  precision  of  the  estimated 
performance  magnitude  j|p  %  as  a  function  of  the  number  of  temporal  samples  M,  the 
number  of  sensor  elements  n  and  a  significance  level  a. 

•  For  N  temporally  correlated  sample  values  on  a  finite  sample  interval  [(7,7' <  <x>]  which  are 

such  that  the  sample  process  covariances  estimates  (Cov{XlXjy,iJ «=  1,2  ,  N  and  i  *  j } 

<,  Nf*'s,  the  derived  form  of  the  Wishart  distribution  G*{||/)-,||,}  -»  Gdl/r’yp-’,  M,  n ) 
asymptotically,  even  for  processes  which  do  not  strictly  satisfy  the  initial  Gaussian 
assumptions. 


REFERENCES 


1.  ESL  Interim  Report,  Kandle,  D.  A.  and  McCarty,  R.C.,  Statistical  Analysis  of  Adaptive 
Beamforming  (ONR  Contract  N00014-85-C-0891),  ESL  Inc.,  Sunnyvale,  CA,  September 
1986. 

2.  Wishart,  John,  “The  Generalized  Product  Moment  Distribution  in  Samples  from  a  Normal 
Multivariate  Population”  Biometrica,  vol.  20A,  1928. 

3.  Wishart,  John  and  Bartlett,  M.  S.,  “The  Distribution  of  Second  Order  Moment  Statistics 
in  a  Normal  System,”  Proceeding  Cambridge  Philosophical  Society  vol.  28,  pp  455-459, 
1932. 

4.  Wishart,  John  and  Bartlett,  M.  S.,  “The  Generalized  Product  Moment  Distribution  in  a 
Normal  System,"  Proceedings  Cambridge  Philosophical  Society,  vol.  29,  pp  260-270,  1933. 

5.  Ingram,  A.  E.,  “An  Integral  Which  Occurs  in  Statistics,"  Proceeding*  ^...mbridge  Philo¬ 
sophical  Society,  vol.  29  pp  271-276,  1933. 

6.  Hsu,  P.  L.,  "A  New  Proof  on  the  Joint  Product  Moment  Distribution,"  Proceedings 
Cambridge  Philosophical  Society,  vol.  35,  pp  336-338,  1939. 

7.  Cramer,  11.,  Mathematical  Methods  of  Statistics,  Princeton  University  Press,  Princeton,  New 
Jersey  1951. 

8.  Fisz,  M.,  Probability  Theory  and  Mathematical  Statists  (3rd  Ed.),  John  Wiley  &  Sons,  New 
York,  1963. 

9.  Anderson,  T.  W.,  An  Introduction  to  Multivariate  Statistical  Analysis  (2nd  Ed.),  John  Wiley 
&  Sons,  New  York,  1984. 

10.  Horn,  R.  A.  and  Johnson,  C.  A.  Matrix  Analysis,  Cambridge  University  Press,  New  York, 
1985. 

11.  Birnbaum,  Z.  W.,  Introduction  to  Probability  and  Mathematical  Statistics,  Harper  and 
Brothers,  New  York,  1962. 

12.  Hildebrand,  F.  B.,  Introduction  to  Numerical  Analysis,  McGraw-Hill,  New  York,  1956. 

13.  Lanezos,  C,  Introduction  to  Numerical  Analysis,  Prentice  Hall,  Englewood  Cliffs,  New 
Jersey,  1956. 


14.  McCarty,  R.  C.  and  Jolley,  D.  K.  Tutorial  On  Prony's  Method,  ESL  TM-1989,  Unclassified, 
ESL  Incorporated,  Sunnyvale,  CA,  May  1986. 

15.  McCarty,  R.  C.,  “Sample  Estimators  for  Entropic  Measures  of  Mutual  Information,"  SP1E 
Proceedings,  San  Diego,  CA,  August  1988.  (to  be  published  as  above) 


ADDITIONAL  REFERENCES.  (NOT  CITED  IN  TEXT) 


Bicnvenu,  G.,  Kopp,  L.,  "High-Resolution  Array  Processing  Using  Eigensystem:  Esti¬ 
mation  Aspect,”  High-Resolution  Spatial  Processing  in  Underwater  Acoustics,  R.  A. 
WagstafT  and  A.  B.  Baggeroer  (eds.)  1984. 

Bienvenu,  G.,  Kopp,  L.,  "Principal  de  la  Goniometrie  Passive  Adaptive,”  Proc.  7'ieme 
Collopue  GRESTI,  Nice  France,  1979,  pp.  106/1-106/10. 

Bienvenu,  G.,  Kopp,  L.,  "Adaptivity  to  Background  Noise  Spatial  Coherence  for  High- 
Resolution  Passive  Methods,”  Proc.  IEEE  Int.  Conf.  on  Acoust.  Speech  and  Signal 
Processing,  Denver,  CO,  1980,  pp.  307-310. 

Bienvenu,  G.,  Kopp,  L.,  “Source  Power  Estimation  Method  Associated  With  High-Rcso- 
lution  Bearing  Estimator,”  Proc.  IEEE  Int.  Conf.  on  Acoust.  Speech  and  Signal  Processing, 
Atlanta,  GA,  1981,  pp.  302-305. 

Cron,  B.  F.,  Sherman,  C.  H.,  “Spatial  Correlation  Functions  for  Various  Noise  Mujsb  ” 

J.  Acoust.  Soc.  Am.,  vol.  11,  1962,  pp.  1732-1736. 

Golub,  G.  H.,  Van  Loan,  C.  V.,  Matrix  Computations,  The  Johns  Hopkins  University  Press, 
Baltimore,  Maryland,  1983. 

Liggct,  W.  S.,  “Passive  Sonar  Fitting  Models  to  Multiple  Time  Series,"  NATO  on 
Signal  Processing,  ed.  G.  W.  R.  Griffiths  et  al.,  New  York  Academic  Press,  1973,  pp. 
327-345. 

Owsley,  N.  L.,  "Sonar  Array  Processing,"  Array  Signal  Processing,  S.  Haykin,  cd.,  1985. 

Stewart,  G.  W.,  "Error  and  Perturbation  Bounds  for  Subspaccs  Associated  With  Certain 
Eigenvalue  Problems,"  SIAM  Review,  vol.  15,  1973. 

Wax,  M.,  Shan.  T-J.,  Kailath,  T.,  "Spatio-temporal  Spectra!  Analysis  by  Eigenstructurc 
Methods,"  IEEE  Trans,  on  Acoustics.  Speech  and  Signal  Processing,  vol.  32,  pp.  817-827 
and  727-764. 

Horibc,  Y.,  “Entropy  and  Correlation,"  IEEE  Trans.  Systems,  Man  and  Cybernetics,  vol. 
SMC- 15,  No.  Sept/Oct  1985. 

Blachman,  N.,  Noise  and  Its  Effect  On  Communication,  McGraw-Hill,  New  York. 

Blum,  J.,  and  Rosenblatt,  M.  Probability  and  Statistics,  Chapter  5,  W.  B.  Sanders,  Phila¬ 
delphia,  PA  1972. 


ESE-88-0099 


Parzen,  E.,  Stochastic  Processes,  Chapter  1,  Holden  Day  Inc.,  San  Francisco,  CA  1962. 

Feller,  W.,  An  Introduction  to  Probability  Theory  and  Its  Applications,  vol.  2,  John  Wiley 
&  Sons,  New  York  1966. 

Grenander,  U.  and  Rosenblatt,  M.,  Statistical  Analysis  of  Stationary  Time  Series,  John 
Wiley  Si  Sons,  New  York  1957 

Parzen,  E.,  Modern  Probability  Theory  and  Its  Applications  ,  John  Wiley  &  Sons,  New  York 
1960. 

Parzen,  E.,  Stochastic  Processes,  Holden-Day,  San  Francisco,  CA  1965. 

Monzingo,  R.  A.  and  Miller,  T.  W.,  Introduction  to  Adaptive  Arrays,  John  Wiley  &  Sons, 
New  York  1980. 

McCarty,  R.  C,  A  Sample  Estimator  of  Beam-Forming  Performance  and  Its  Distribution 
Function,  ONR  Report  ESI. -87-0317,  ESL  Inc.,  Sunnyvale,  CA  1988. 


11-2 

''-v  +-  *  -*  r  /  j-  f  /  a  f  -j-  V  Vj.  -j.  -j.  v  V  V  'J  ’j>  "J "J  *-•  V  '  ■  - 


w  \« ■> -  -  •  ■> V 


I 


» 


> 


I 


§ 


I 

.J 


S 

$ 

I 


